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Abstract 

We establish a framework to construct a global solution in the space of finite 
energy to a general form of the Landau-Lifshitz-Gilbert equation in R^. Our char- 
acterization yields a partially regular solution, smooth away from a 2-dimensional 
locally finite Hausdorff measure set. This construction relies on approximation 
by discretization, using the special geometry to express an equivalent system 
whose highest order terms are linear and the translation of the machinery of lin- 
ear estimates on the fundamental solution from the continuous setting into the 
discrete setting. This method is quite general and accommodates more general 
geometries involving targets that are compact smooth hypersurfaces. 

1 Introduction 

Micromagnetics, a model based on the work of Landau and Lifshitz, offers a description 
of magnetic behavior in ferromagnetic materials. This continuum model is based on 
finding critical points of the micromagnetic energy associated to the magnetic moment 
represented by the field u : fl —>■ S^: 

E[u]= / |Vm|^ + k / / \VM\'^dx-2 / h^t ■ u 

Jn Jn Jr^ Jn 

where Q is the region occupied by the ferromagnet and all the physical constants have 
been normalized to be 1. These four terms are known as the exchange, anisotropic, 
magnetostatic, and external energies respectively. The exchange energy penalizes spa- 
tial variation over macroscopic distances. The anisotropic term favors special directions 
of the magnetization. The magnetostatic energy represents the work required to build 
up the macroscopic body by bringing its magnetic moments from infinity to their even- 
tual position; M is nonlocal and is defined by AM = div u where u is extended to be 
outside of Q and the equation is to be understood in the sense of distributions. Finally, 
the external energy favors magnetization aligned with an external applied field. The 
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first variation of the energy witli respect to u is called the effective field: 



5E[u] 



6u 



Constrained to the set of candidate solutions {u : Q S^}, the resulting dynamic 
equations associated to this energy are, in dimensionless form, given by 



where A denotes the wedge product in M^. To incorporate the Gilbert damping law, 
whose origin lies in the observation that such systems reach equilibrium and must 
have decreasing energy over time, a dissipative term can be added on resulting in the 
following equation: 



In the physical setting, the magnetization configuration is the compromise that results 
from the competition of all of these terms to minimize the energy. However, many 
relevant features of solutions to the full equation are captured by retaining the term 
associated to the exchange energy in which case /ig// = ^u. For the question we 
address here, in particular, the other terms can be considered lower order terms. An 
analytically interesting variant of the problem results from attempting to enlarge the 
class of target spaces beyond S"^ to hypersurfaces with unit normal u. Although this 
variant does not hold common currency in micromagnetics, this might be a reasonable 
means of incorporating anisotropic contributions. 

The Cauchy problem for this generalized Landau-Lifshitz-Gilbert equation (LLC) 
is the problem of finding u, given initial data / : ^ A/" for A/" a compact smooth 
hypersurface in M.^ with normal vector u, satisfying 



where a > 0. We will refer to the first term as the Schrodinger term (this terminology 
will be clarified later) and the second as the damping term. When only the damping 
term is present, this equation is the harmonic map heat flow problem, which will later 
be a focus of discussion. Although standard usage of LLG refers to in the special 
case A/" = 5"^, we adopt here a usage to accommodate this general class of target spaces. 

There are several standard forms of which are equivalent for smooth solutions. 
The first results from the vector identity A {v AC,) = — f + {v,C,)v which holds for v 
a unit vector. From this, (jl.lj) can be written as 



dtU = uA heff, 



dtu = u A hcff — au A {u A /ic//)- 



dtu = u{u) A Am - au{u) A [uiu) A Au) 
m(x,0) = / 



(1.1) 



\u = u{u) AAu + a{Au + (Vz/(m) 



Vu)u{u)) 



(1.2) 



From (jl.l|) and (jl.2p . two additional formulations follow: 



adtu - u{u) A dtU = (1 + a^){Au + (Vi^(n) • Vu)ij{u)) 



(1.3) 
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dtu + ai'{u) A dtu = (1 + a^)i'{u) A Am 



(1.4) 



The question that we address in this paper is whether global solutions exist in the 
class of finite Dirichlet energy. The main result in this paper is the following: 

Theorem 1.1. For any f G H'^ {M."^ , Af) , there exists, for t > 0, a solution to / ti.i)) 
which is smooth away from a singular set that has locally finite 2-dimensonal HausdorfJ 
measure with respect to the parabolic metric. 

Weak solutions, even partially regular ones, have been shown to exist for such 
equations in the case Af = S"^. The techniques developed, however, often exploit the 
special geometry of S"^ and side-step the natural difficulties that the general problem 
poses. Consequently, these techniques are not easily adaptable to the setting that we 
consider here. One category of results showing existence of weak solutions has made 
crucial use of the fact that LLG can be written in divergence form in the case A/" = 5*^; 
this in turn gives rise to a definition of weak solution. Amongst these is the work 
of Alouges and Soyeur [Jj who have shown that in three dimensions, energy bounds 
are sufficient for existence of such a weak solution. Guo and Hong [3^ successfully 
carried through the argument that Struwe in [7j employed for the harmonic map heat 
fiow to exhibit a Struwe solution, a partially regular solution that satisfies an energy 
inequality and is smooth away from a finite set of point singularities. In this dimension, 
Struwe shows that the existence of smooth solutions in addition to energy estimates 
is sufficient for the characterization of such a partially regular solution. In the case 
of the sphere, the problematic term u A Au does not preclude energy estimates since 
multiplication by u and Am eliminates this term altogether. In addition, they show 
that the linear term in LLG is strictly parabolic and have used parabolic theory 
to establish local existence. It is worthwhile to note that uniqueness in the class of 
partially regular solutions is open; our partially regular solution may be different from 
the Struwe solution. Through private communication, we have learned that Melcher in 
[Uj has recently attained a characterization of partially regular solution similar to ours 
in three dimensions for A/" = S*^ using a Ginzburg-Landau approximation. 

In the context of existing results, the main difficulties are already apparent. The 
nonlinearity in the highest order z/(m) A Au term poses the most apparent obstacle. At 
first glance, it renders unavailable the machinery of partial regularity theory established 
for semilinear equations: linear estimates on the fundamental solution (e.g., — L'^, 
Strichartz-type) and regularity arguments based on the inversion of the fundamental 
solution. An additional problem that this term poses for general target Af is that the 
local existence of smooth solutions is not just an immediate consequence of parabolic 
theory. Crucial tools that are lost include the maximum principle as well as a Bochner 
identity. Lastly, the divergence structure of the equation is lost when we depart from 
the sphere. 

Our method is based on approximating systems on a uniform spatial grid. The idea 
is to show that a sequence of solutions to this approximation converges to a solution 
to the original equation as the grid-size approaches zero. The idea of approximating 
systems via discretization is not new; in fact, one of the earliest works investigating 
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for the case a = 0, Af = S'^ jH] uses discrete approximations to show local exis- 
tence. However, extendings its scope to partial regularity is new. While the techniques 
developed here are demonstrated for the specific setting of the LLG equations, the 
generality of these techniques should be emphasized. This key ingredients used are: (i) 
the construction of a suitable discretization of the system; (ii) establishing a suitable 
linearization of the discrete system; (iii) the construction of the discrete fundamental 
solution of the resulting linear operator and deriving appropriate linear mixed space- 
time estimates. The setting that enables the discrete constructions described by (i) 
and (iii) are the subject of Section |21 

The advantages of this method address the difficulties of the problem. By dis- 
cretizing the system, existence of solutions to the resulting ODE is immediate once 
we attain uniform bounds; there is no need to appeal to heavy machinery to conclude 
that smooth solutions exist. The transformation which results in a linear highest or- 
der term opens the door to the arsenal of linear estimates and methods. The explicit 
construction of the discrete fundamental solution and the resulting linear estimates 
permit higher derivative bounds that has no dependence on the special structure of 
the sphere. 

Before employing the proposed method to construct a partially regular solutions for 
LLG, we illustrate in Section El this method for the harmonic map heat flow problem, 
which involves only the damping term. The same characterization of a partially regular 
solution has been achieved by other methods 2, , but we consider this excursion to 
be worthwhile since it is a setting in which we can isolate the difficulties that just one 
term poses (makes use of ingredients (i) and (iii) above) before adding on the separate 
difficulty involving the Schrodinger term. In Section EJ our main result is proved. 



2 The Discrete Setting 
2.1 Notation and Definitions 

Consider a uniform grid on M'^ determined by the lines Xi = jih, i = 1,2, ... ,d where 
ji are integers. For vector valued functions defined on the grid, the definition of 
scalar product and discrete norms that will be used are: 

3 

\\A\Li = {h'Y.\^)\n"'. 
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where = = u{jih, . . . ,jdh) and = u{ji, . . . , jih + h, . . . ,jdh). The basic 

difference operations are 



and 

The discrete Laplacian A'' can be defined as 

i=l 

For some multi-index a — {ai, a^, . . . , a^, a^), ctj^, e N, we will use the notation 

d 

i=l 

and 

i=l 

The discrete Sobolev spaces are defined as 

|a|<it 

with = W^'"^. Their homogeneous counterparts are defined by 

|a|=ifc 

with Hh'' = Wh'"^. 

For integer k, we will use the notation D^Uj to denote the vector consisting of all 
difference derivatives D"'Uj, \a\ = k. For the case k = 1, and D~ will denote the 
approximate gradients, the vectors with components D+ku'j and D^^Uj respectively. 

Notable Imperfections of Discretization 

Discretization introduces some imperfection which are important to highlight. The 
common operations of product differentiation and integration by parts have a slight 
shift. The forms of product differentiation that will be used are: 
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For functions u^, both belonging to L^, the summation by parts formula- the discrete 
analog of integration by parts- is given by 

Another important imperfection that arises is that difference derivatives lose their 
'tangency'. Consider a map u from to A/", where A/" is a hypersurface in with 
unit normal v. Difference derivatives of u no longer lie in the tangent plane Tu-M. 
To negotiate some of the difficulty that this defect creates, we will often separate the 
portion of the difference derivative that lies in T^-M and that which is orthogonal. 

2.2 Basic Inequalities and Some Useful Tools 

Most of the inequalities that will be used have long been established. Amongst these are 
Holder's, Minkowski's and the U' interpolation inequalities which translate perfectly 
to the discrete setting. Discrete versions of the Sobolev embedding theorems were 
established by Ladyzhenskaya jH] using interpolation operators. In particular, we will 
make use of the following discrete Sobolev-interpolation inequality: 

\\D^u%Lr < CWuXAlD^u^^, --L = + e{- - -). (2.1) 

The interpolation operators that supply the proof for the above inequalities are useful 
and we mention one which will play an immediate and a later role. Confining our 
attention to the case d = 2 for the remainder of this section, associate with u (which is 
only defined on the uniform grid in the plane) the unique interpolated polynomial ph 
which matches the value of u at each grid point and is of the form ao + aix + a2y + asxy 
in each square. Specifically, for j indexing pairs (ji, ^2), 

Ph = ^iuj + D+iUj{x-jih) + D+2Uj{y-j2h) + D+iD+2Uj{x-jih){y-j2h))xn,ix,y), 
j 

where Xaj{x,y) denotes the characteristic function on the square Kj. Across any one 
side of a square, ph is linear, so we have that ph is continuous across interfaces, ph is 
an iJ^-interpolant associated to u. Of particular importance are the equivalence of the 
norms ||ph||Lp and Hm'^Ulp as well as || VphlUp and for any p > 1. Equipped 

with this interpolant, we can prove the following propositions. The first is a localized 
Sobolev-interpolation inequality. 

Proposition 2.1. Let ^] C and ( e C^{^). Then 

Remark 2.2. Without the (, the quantity of interest is \\u^\\'^]^4b(^2-^ o,nd the proposed 
inequality is given by 1^2.1^ for n = 2,r = 4s, j = 0, ^ = 1,6* = ^. 
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Proof. This inequality follows from the equivalence of the norms and Hi^'^IIlp 

as well as ||Vp?i||lp and ||D^m'^||/^p for any p > 1 and the equivalent statement in the 
continuous setting, as proved in Appendix A. □ 

Proposition 2.3. // ||m'^||j:^2(-jj2) < C independent of h, then there is a subsequence 

{P/ifc} that converges strongly in H^{M?). 

Proof. Notice that the boundedness of ph in immediately gives weak convergence 
in H^. To show strong convergence, we employ Ascoli-Arzela. In particular, we 
must show that {Vph} are bounded and equicontinuous in LP'. 

We have LP bounds on Vpu by construction. To show equicontinuity, we must show 
that Ve > 0, 35 = ^2) = (/ci/i, k2h) such that \\Vph{- + (^1, • + ^2) - Vp/i||i2 < c\5\. 
We have 

||Vp/.(- + 5i,- + 52) - Vpfe||i2 

1+fel J2 

D+iD+2Uj^,j^){y - j2h)\'^ + |(D+2tiii,j2+fe2 - D+2Uj,,j^) + {D+iD+2Uj,,j^+k^ - 
- D+iD+2Uji,j2){x - jih)\^dxdy 

< ^ X/ / l-^+i'"ii+fcij2 - -0+i'"iij2l^ + l-D+2«iij2+*;2 - -D+2'"ii,i2r + 
h,32 ''^n'i2 

\D+iD+2Uj^+kuj2 - D+iD+2Uj^,j^{y - 32h)\^ + \D+iD+2Uj^,j^+k^ - 
D+iD+2Uj,,j^{x - jih)\^dxdy 
= 2{L + LL). 

To bound /, we can write the difference D^iUj-^+ki,j2 ~ -0+i^ii,i2 

ki 

^+l'^jl+kl,j2 ~ -^+1%1,J2 — ^(-^+I%l+fej2 ~ -^+l%l+fe-l,J2)- 

fc=l 

Since (oi H h a„)^ < 2n(a^ H h a^), we have the bound 

ki 

\^+i'^ji+ki,j2 ~ D+i'^ji,j2f ^ 2/ci ^ \D+iUj^+kj^ — D^iUj^+k-i,j2\'^- 

k=l 

Similarly, we have the bound 

k2 

\D+2Ujij2+k2 - -0+i%i,j2l^ < 2A;2^ \D+iUjj^j^+k - D+iUj^^j^+k-i]'^- 

k=l 

Using the given information that ||w'*||ij2 < C, we have in particular that 

h 

W^^Wk = E \^+^^hJ2 - D-,Uj,jf + |i5+2<,,2 - D-2uljf < C, 
ii,i2 



which allows us to conclude that 
I = 



/ l-O+l^i+feiJ^ - -D+l^jlJ2r + |-^+2'"ji,j2+fc2 - D+2Uji, 

jl,j2 "^°ilJ2 

fcl 



jlJ2 fe=l 

2(^2-1) ,j2+fc-l| J 



jl,j2 k=l 

,2r 



< 2h'[{ki - l)kiC + (k2 - l)k2C] 

< 2Ch\kl + kl) 

< c\5\\ 

By assumption, we also have that "^j^j^ ^^l-0+i-02%ij2p < C* so 

|-D+l-D+2liji,i2+fc2 ~ -0+1-0+21*^1^2^ 
< Ci/l^(fci(fci- 1)^ + ^2(^2-1)^ 



□ 



2.3 Constructing Discrete Approximations 

We now construct semi-discrete approximations for LLC The method that we adopt 
is that of finite-differencing on a uniform spatial grid with gridsizc h. This method 
reduces the problem to a system of o.d.c's with the unknowns being {m^}, the value of 
the function at every grid point. We seek a construction that yields a stable scheme. 
A natural way to achieve this is to start by discretizing the energy rather than by 
discretizing the equation. Exphcitly, the Dirichlet energy can be discretized to yield 
the discrete energy E'^[%jt\: 

j 

where j indexes over the grid pairs {31,32)- The discrete analog of is 

^ <+i - 2^? + <i ^ _..n 
5u) h? ^ 
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The discretization of the system associated to this energy is then given by 

\ u){t = 0) = f^ ^'-'^ 

where denotes the unit normal vector to A/" at Uj{t) and the Lagrange's multipher 
X'^{t) is given by X'j{t) = -A'^u){t) ■ Vj{t). 

We also have the discrete equivalent of ()1.2|) : 

adtu%t) - Vj A dtu] = (1 + a2)(A\^ + \)vj). (2.3) 

This discretization yields energy bounds. In particular, since Vj is orthogonal to 
we can multiply ()2.3|) by and sum across indices j to get 

d,E\u\t)\ = {d,u\ = -c^WdtuXl (2-4) 

Since a > 0, we can integrate on any time interval [0, T] to get the discrete energy 
bound 

E^[u^[T]] < E^[f ]. 



2.4 Constructing Discrete Fundamental Solutions 

Our method of proof is based on linear estimates. Although not immediately apparent, 
the relevant linear equation to consider is the linear damped Schodinger equation. Since 
we will also be looking at the harmonic map heat flow, we will also consider just the 
linear heat equation. We now construct fundamental solutions for the discrete analogs 
of these equations. The fundamental solution of the discrete heat equation is given by 
$^(t) = where solves the system 

$?(t = 0) = 

The notation is consistent with the definition of the weighted h norms that is being 
used and yields || )^5''||L2(Md) = 1- An explicit description of the fundamental solution 
follows from the observation that the discrete Laplacian A^u'i = A^u^ , . can be 
written as the sum of operators. 



=1 



As in the continuous case, the description of the fundamental solution in higher di- 
mension rests on the one dimensional problem of solving 



^n^t = ^) = ¥1 



(2.5) 
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The kernel associated to the operator ^ — j^D~^ is given by K'^, where 

otherwise 
and that associated with ^ + j^D~ is given by K'^ where 

[ otherwise 

we exploit once again the commutativity of and to conclude that the solution 
to ()2.5j) is given by the discrete convolution 

The fundamental solution of the discrete heat operator in M*^, $j , is then just the 
product of the w^.'s over 1 < i < d. 

Similarly, the fundamental solution of the discrete damped Schrodinger equation, 
denoted by W"''* = (U'^'^)j, solves the system 



U^'\t = 0) = 



The construction proceeds exactly as in the case of the heat equation. The following 
proposition summarizes these constructions: 

Proposition 2.4. 

1. = the fundamental solution for the d- dimensional discrete heat equa- 
tion, is given by 

1 { _L\k-ji { t \k 

The solution operator to the linear discrete heat equation for inital value f^ is 
given by the discrete convolution <l>'*(t) * f^ and will be denoted by (^^{t)f^. 

2. 14°"'^ = (U'^'^)j, the fundamental solution for the d-dimensional discrete damped 
Schrodinger equation, is given by 



The solution operator to the linear discrete damped Schrodinger equation for inital 
value f^ is denoted byU°''^{t)f^. Additionally, denote U^'^ = U^. 
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2.5 Discrete Linear Estimates 

Having constructed the discrete fundamental solutions, it is possible to derive discrete 
analogs of the well known U' — L'^ estimates. 

Proposition 2.5. (L\ — L\ Estimates) Let KJ^ = $^ orUa,h, a > from Proposition 
Then the following estimates hold: 



ll/^ll-L' 

\f^^it)f''\\Ll{R'i) < C d (2.9) 



(1) 



l^'/C^W/'lUsfM^) < c (2.10) 



(2) 

\\D']C.'^(t)fHr.,„..<C- 

Proof. We will show the derivation of both estimates by taking /C^ = Since a > 0, 
the same derivation holds in the case that JC^ = lA"'^. 

These estimates will follow upon obtaining suitable L\ and estimates on $^ 
and Z)^$^. Specifically, if the bounds 

ii$'^(t)iu. <i, \mt)\\,^<\ 

t2 

hold, then by Hausdorff- Young we have 

II*'(^)/'IIls < ll*'(t)IUdl/'llL?, 



where 1 + i = i + 1. By interpolation, ||$'^(t)|Uj; < {\\<^>\t)\\Lif-%\\<^\t)\\L^y for 
= 1 - i = 1 - (1 + i - i) = i - i, so \\<^\t)\\Lr^ < from which (1) follows. 

Similarly, (2) follows from the estimates 



1 I d 



t2 t2 + 2 

To show ||$^(t)||Li < 1, it suffices to get an estimate on the kernel of the one 
dimensional problem, Wj-. (t) given by (|2.6|) which is a discrete convolution of two 
functions in L\. Since the L\ bound of a convolution of two L\ functions can be 
bounded by the product of the two L\ bounds, we have 



^ 1 



'_L)k-ji (^)k 



_ 2t 2t 

< 1. 



Since w^^{t) is a convolution of two L\ functions (in jj) , the bound of w'j.it) 
will be the product of the bound of any one of the convolved pair with the 
L\ bound of the other. Using Stirling's formula, and re-indexing with k, we have 



Writing x = -A, and letting f{k) = (f )'^e we have 



—L — ( ( — \^\' p~^+^ _|_ (^\k -x+k 
dk ' 

k k k 
= e-^+^f-)Mn^ 



Thus / is maximum when |^ equals which occurs when fc = x= p-^^ = -y/^. 
Returning to the bound of w^(t), we have 

\\wl{t)\\,^ < CW^fmL- < C||-^|U^/(x) = ci= (2.11) 
as needed. It follows that ||$''(t)||L-' < 

To get gradient bounds on the fundamental solution, we express 

For bounds on the gradient, we use the fact that the gradient will only 'hit' 
one term in the convolved pair and that the other is in L\. Again, let x = -p, 
and f{k) = [^)^e~^~^^. Since f{k) achieves its maximum at k = x, 

Using the bound given by ()2.11|) for each of the other d — 1 Wj_ (t) terms in 
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the product, we get 



2 t 2 + 2 

Therefore, 



C 

< - 



t2 + 2 

Employing again the fact that the L\ bound of a convolution can be bound 
by the product of each of the L\ bounds in the convolved pair (and hence 
||D+aW^^(t)||^i < ^ < for each a), and that each of the Wj.{t) is an L\ 
function , we conclude that 

< C sup \\{D+,w^^{t))Ui^,w^^{t)\\L^ < ^. 

l<i<d t'i 

□ 



3 Harmonic Map Heat Flow 

The first problem on which we apply the machinery that has been developed so far is 
the Cauchy problem for the harmonic map heat flow problem in two dimensions. This 
is the problem of finding m, given / : — A/" such that 

dtu = Au + (yu{u) ■ Vu)v{u) , . 

uiO) = f ■ 

Af, as usual, is taken to be a compact smooth hypersurface. The semi-discrete ap- 
proximation for this equations follows the construction in Section 12.31 and is given by 

where I'j = i^{u^) and the Lagrange's multiplier A^(t) is given by 

\%t) = -A^u){t).v,{t). 

In what follows, the letter C will designate a generic constant that may depend on 
M but is independent of the gridsize h or any solution unless specified explicity like 

Estimates on 

For the purpose of attaining bounds on solving ()3.2p . it will be useful to preface 
the technical detail with several observations concerning A^. In the continuous setting, 
there is an explicit expression for the Lagrange's multiplier due to the fact that diU G 
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TuM . This permits the following manipulation: 

—Am ■ ^{u) = di{—diU ■ I'iu)) + {diU ■ div{u)) = Vz/(m) ■ Vu. 

Thus, 

|Aj| = \ iyu{u) ■Wu)\< C|Vn|l (3.3) 

In the discrete setting, difference derivatives lose their tangency which robs of as 
clean an expression. Writing a) = D+fcM^ ■ z/,-, we can attempt the same manipulation 
as in the continuous case to get the following expression for A^: 



Aj = D^.a] - {D^iU ■ 



The first observation is that even though this expression has a local defect, we can still 
establish a discrete analog of ()3.3j) . 

Observation 1 

|AJ| <C(|Z}+,n,f + |D_,^i,f), (3.4) 

Expressing 

we need only show that for fixed index fc. 



Fix some 6 > 0. Then if ~ '^'jj — '^'] have 

l^'?- — I 

in which case we can take C = ^. For the case l^j^^-^ ^ u'jj < 6, the desired inequahty 
is a direct consequence of the compactness of A/". Namely, since A/" is smooth, we can 
express A/" locally as the graph of some quadratic function, so there is a ball about 
every point such that M behaves quadractically. Due to the compactness of J\f, there 
is a finite subcover of these balls. If we choose 6 to be the minimum of these radii, and 
to specify henceforth that the gridsize h is to be smaller than S, then the interpretation 
of (Mjj.^j — ■ Uj as the projection of — Uj^ onto uj immediately yields the given 
inequality. 

Observation 2 

For 6 > 0, choose h such that for k = 1,2, — < S, — u'j^,-i\ < Then 

\D'\^\ = 0i\D'u1\MD'u1\\D\1\), (3.5) 
where / indexes over next to nearest neighbors of j. 
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This follows from the calculation 

= o(iDVlM^Vll^Vl), 

I indexes over next to nearest neighbors of j. 



3.1 A Global Smooth Solution for Small Energy Data 

The first result that will be shown is that under the assumption that the initial data 
has small total energy, there exists a smooth global solution. While the result in this 
section can be seen as a corollary to the proof of the finite energy data case, this 
assumption permits the 'fastest' path to our goal . The process is laid out here to fix 
ideas and to motivate the line of reasoning underlying some of the extra steps needed 
in the next section. The following result will be proved: 

Theorem 3.1. There is a constant such that for G -ff^(]R^) andsnpf^E^ lf'^] < Cq, 
the associated sequence of H^-interpolants {ph} converges strongly in to p, which 
solves ij) in the sense of distributions. 

Proof. The proof of this result is a two step process of attaining appropriate bounds 
on the sequence and using these bounds to construct interpolants that converge to 
the solution claimed. 



Step 1: Attaining Bounds 

The only quantities on which we have bounds immediately are energy and dtu^. Namely 

E^K(T)] + ||a,^1li.([o,T],L|) < (3-6) 

Define 

= {m'^ = {u';}j\t^/^\\D\\t)\\Li < oo Vt}. 

The higher derivative bounds that we aim for here are uniform- in- /z. bounds in (i) 
L^Hl), (ii) T'^ and (iii) Lip(ifi). 



L (H^) To get bounds on terms involved in ||Z} m 11^2, it suffices to bound 



IIA'^M^II r2 since the discrete Laplacian controls all second difference derivatives. 

h 

Using Observation 1 and (|2.ip for j = 1, r = 4,p = g = 2, and the global energy 
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bound dSini) 



< 


C{ 


\dtu\t) 


r2 

'^h 


+ (/i'S,|Ajn2) 


< 


C{ 






1 1 1 1 7^ 1 h { a\ 1 2 1 1 \ 

+ 111^^^ WI IIl^) 


< 


C{ 


\dtu\t) 






< 


C{ 


\dtu\t) 






< 


C{ 


\dtu\t) 




^ E\f\h\\DW\t)\\^.^ 


< 


C{ 


\dtu\t) 




+ e|pV(t)|L.) 



Now choosing eo sufficiently small, we can absorb the C(eQ )||D^m'^||^2 term into 
the left hand side to get 



Now integrating in time, we can make use of (j3.6|) again to conclude that 



To get bounds involving D'^u^ bounds, we can take difference derivatives of the 
equation and and then gain a derivative by using the estimates on the derivative of 
the fundamental solution. Using the estimate given by ()2.10|) , once with p = q = 2 
and another with p = 2, g = p' = 4/3, 



Observation 1 and 2 on Ay together with the smoothness of u yield the inequality 



D\\)uj) < C\iD'X']\ + \X^\\D'u^\) < C{\DW^ + \DW\DW), 



where / indexes over next to nearest neighbors of j. We can now apply ()2.H) in 
addition to repeated use of the discrete Holder and Minkowski inequalities to get 



(ii) JF^ For each index j, u'j satisfies 
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Proceeding, 

< c(PV'IL|) + ce'' ^^-^{1|dV(s)1L41|dV(s)1l. + WD'u^Wl.^ds 

t>0 Jo (t — s)-'' h 

Letting = supo<s<tS^^^\\D^u''{s)\\L2 



Noting that 



y{t) < e'J' + c{eo)y{trh^" / — ^ 



and 

^1/2 



we have that ?/(t) satisfies 

y{t) < Cef + C{eo)y{tf'^ 

for all t G [0, oo). The function L{y) = y — C{eo)y^^'^ has one root at ?/ = and 
another positive root, reaching its maximum at y = g^^^- Now we can choose 

eo so that Ce^^ < L(g^j^^). Since t i— > y{t) is continuous and y{0) = 0, we can 
conclude that y{t) is bounded, uniform in h. 



m 



Lip{Hl) The energy bound ()3.6j) immediately implies that m'^ G L^{H\). We 
can do better using the bounds that we just attained in T^. For ti < t2, let 
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t* — ti + ^^y^. The bound we seek is a consequence of the following: 

\\D\u\t,) - u\m\Li 



'tl 

rt2 



+ 



1 1 




_5|l/2l 




1 1 




-S|V2I 






\t2 


_5|l/2 


f" 


1 


hi 





\\D'u\s)\^.\\D\\s)\^.ds) 
1 1 1 



< Cieo)iit*-h) + ih-t*)) 
= C{eo){t2-h) 

Step 2: Convergence of the Interpolants 

We have shown that there exists, for each h, sequences of solutions — {uj{t)} such 
that 

with uniform bounds in h: 

sup(||h'^IL.(^.) + supi^/li^V(i)|U. + sup ——- ^) < C{eo) 

h ^ t ^ h^t2 \t2 - tl\ 

Associate with each sequence {u'j{t)} the interpolant Phif)- Then for (f) e C^(R+ x 
M^), we will show the following statement: 

- / Ph-dt(t)+Ph-^(t>+{'^Ph-yy{Ph)HPh)-(t>= P/^(0) • 0(0) + 0(/i) (3.7) 

Recall that, for {x,y) e Dj, j — (^1,^2), the interpolant is given by 

Ph{x,y,t) = u%t)+D+^u]{t){x-j^h)+D+2U%t){y-j2h)+D^^D+2U%t)^^^ 

For j = (ji,j2), the notation /□. will denote We will also slightly 

abuse the notation (which is only defined on grid points) to play the dual role of the 
function which takes the value of in Hj. Since e LP'{Hl)nL°°{H\), the following 
are immediate: 



(i) / (/ p^.dt4>-Y. / «?-^t0) = O(/i) 

(ii) / (/ p,-A0-^ / u) ■ ^ct>) ^ 0{h) 

Jo JR2 ■ JDi 
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From (i), 

POC P POO P P 

Jo J Jd, Jo j Ja, j Ja, 

Using the equation ()3.2|) for dfUj, what remains to be shown is that 

j 3 
POO P 

(i^) / E/ A>,-0+(Vp,-Vz/(Ph))KPh)-0 = O(/i) 

^0 , JOi 



J 

To show (iii), it suffices to show that 



JO JD, ^0 , Jn, 



^. 



' h 



, Ja 



/•CXD /• 



What remains is (iv). Associate with the step function = (pjXOj- We treat this 
term by sphtting the integral into two parts: 







and 



POO P 

= / E + (vp^ ■ vz/(p,))z/(p,)) ■ ■ 



Since G L^(if^), / = 0{h). The reason behind the sphtting is really to negotiate 
the local defect regarding A^. Fixing j, \^ does not come close to cancelling out the 
relevant terms in Vph{jh) ■ Vuijph^jh)). However, just as it was seen in Observation 
1 that the local defect did not prevent us from attaining a global statement that was 
satisfactorily close to the continuous case, we can also show that by considering the 
sum of the two terms in // as two sums, that the terms that result have a satisfactory 
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cancellation. The only thing that was preventing us before from switching the spatial 
sum and integral was 0. Replacing by resolves this immediate problem. Keeping 
in mind that the terms in Vph{jh) ■ Vz/(p/i(j/i)) that require cancellation are of the 
form D+iu'j ■ we consider the sum J2j ^^^j ' 0^- Letting a* = {D+kU^ ■ Uj), 



j j 

= - ■ D+iiyj)iyj^+i ■ 0j+i 

j 

This term is already in a form that the cancellation is apparent; the local defect in Xj 
has essentially been transferred to a local defect in i> and 0, but these we can handle 
because they are smooth. 

rcx> I* 

= V0,-{Ajz/,+D+,wJ)K),z/(p,)} + O(/i) 

Jo JUj ■ 

POO 

Jo ^ 

POO 

= 0{h). 

(i) - (iv) together entail (j3.7|) . With this expression, we can split the time interval 

into subintervals [0, 5) and [5, oo). On [5, oo), we have uniform if^ which implies from 


Proposition 12.31 that there exists a subsequence {hk} such that Ph„it) p{t). For any 
hk, then, we can use in addtion (i) - (iv) and the fact that G Lip(Hjj) to show that 



Jr^ 



Ph^ ■ dt<P + Ph^ ■ A0 + iyph^ ■ Vi^{ph,)MPhJ ■ 

ph^: POO p 

= (/ +/ ) Vh,-dt(t) + Ph,-^(t>+{Vph,-Vv{pj,^))v{pj,^)-ct) 
Jo Jhk Jm? 

POO P 

= Cieo)hk + Oihk) + / Ph,- dt<p + Ph^ ■ A0 + {Vph^ ■ \/u{phJ)u{ph,) ■ 

POO p 

/ / p • 9t0 + p ■ A0 + (Vp ■ Viy{p))iy{p) ■ 
Jo Jm.^ 

from which we can conclude that p solves ()3.1|) in the sense of distributions. □ 
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3.2 A Partially Regular Solution for Finite Energy Data 

Now we consider the discrete system only assuming data. While we do not have 
total small energy, our initial data is still of finite energy. Locally, then, the energy of 
the initial data can be made small. What is needed is a localized version of the "small 
energy =^ regularity" statement that was gotten in the last section. Achieving such 
a statement will allow us to estimate the size of the concentration sets which in turn 
will allow us to extend the associated interpolants across these sets in such a way that 
they will converge as /i — in the sense of distributions to a solution of the continuous 
problem. We proceed now to show the following result: 

Theorem 3.2. For any uq e H\R'^,M), there exists a global, partially regular solution 
to smooth away from a locally finite 2-dimensonal Hausdorff measure set with 

respect to the parabolic metric. 

The core statement to be shown is: 

Lemma 3.3. There is an cq > such that if E^[u^(to); B-j^^xq)] < cq then 36, inde- 
pendent of h such that Vr, to < t < + ^R^, 

sup ||DV(t)lU2(^_(,„))<C 

h -g- 

T<t<to+STf 

Notation : E^[u^{t); Bfi{xo)] denotes the energy ofu^{t) indexing over j = (ji,j2) such 
that jh E Br{xo) . 

Proof. The goal is to achieve uniform bounds in a reduced parabolic cylinder. We 
will pass through the following intermediate stages to achieve this end: 

(i) Show a local energy inequality 

(ii) Show G L'^{H^) in a reduced parabolic cylinder 

(iii) Show D^u^ G L'^lLf^) for r < oo in a reduced parabolic cylinder 

(iv) Show G L'^{Hl) in a reduced parabolic cylinder 

Proceeding, 

(i) For any e > 0, E''[u''{to); 52/j(xo)] < e/2, 36 such that Vt, to < t < to + 5R^, 

snpEVity,Bnixo)] + lldtuXnitoMLliB^M)) < (3-8) 

h 

Pf: Denote by C, the cut-off function that equals 1 for jh E Br^xq) and for 
indices outside of B2r{xq) with \D^C\ < "I- Multiplying (j3.2p by and 
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summing over j, 

j 

j 



du''- 



'L'i 



h 



. ^ + l^Y. 4.^1: 

3 j j 

Since the same expression holds with D+j replacing D_i, we make use of the total 
energy inequality to get 

Integrating from to to t, 

In particular, choosing 6 = 2E^[fh] i S®^ desired result. 



(ii) L'^{H^) For e > 0, 35 such that 

E^[u\t)- B2r{xo)] < 6 ^ sup < C{e) (3.9) 

P/; To get bounds on terms involved in ||D^ii'^||^2 , it suffices to bound || A'*k'*||^2 
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since the discrete Laplacian controls all second difference derivatives. For any 
cutoff function ( and solving (|3.2|) . UjQ satisfies the following equation: 



A'iu^Q) (3.10) 
-1 
4' 



+ ^K+i + 2«J. + <-i)A%} (3.11) 

= {^(0.+i + 2a + 0.-i)(^<-A?^,) + 

+ 2(D+,0 + /^-.0)(^+.4 + ^-O + i(4.+i + H + <-i)A'C} 
= I + II + III + IV 

This is just the discrete analog of the expression 

A(mC) = AuC + VmVC + uAC = {dtu - Vu ■ Vz/(m))C + 2VmVC + uAC. 

Choose ( such that it equals 1 at indices j such that jh G Br^xq) and outside 
of S2h(xo) with \DW < I and \ A''C\ < Then 

- \\ii\\Li^^.)<c\\\D\\t)\H\\Lim 

- WlllUum < C\\\D\\\D\\t)\h2^^.) < ^\\D\\t)\\,.^s,ni.o) < ^ 

- l|/^llLj^{M^) < C|||A'^Cll«'(t)lllL|(R2) < ^ 

where III is bounded by using the local energy inequality that we attained in 
(i). Now integrating on an interval of length 6R^, we can see that III and IV 
are lower order terms in the sense that 

~ \\m\\L^{[toA-,+5R'2]{Ll{K2)) < C{t)5R 

- \\iy\\L\[to,to+SR'^](Llim) -^^ 

We plan to use the local space-time bound on dtu^ from (i) to handle /. II, 
however, poses an obstacle since there is a 'missing' ( term which prevents the 
step that we used in the small energy case: 

II I ni^>|2|| ^ II ni„>i|2 
\\\U U r2 < \\U U ta 
II I I ll-t.^ — II ll_L^ 

The extra work that is needed is performed in Proposition 12.11 corresponding to 
s = l,p = g = 2. Getting estimates on D'^{u^(), and denoting the lower order 
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terms, l.o.t. = ||/J/||i2(K2) + 11^2(^2), we have 

< C||A(n^t)C)|U2(M2) 

< C\\\D\\t)\'C\\L.^^.) + \\dtu\t)C\\L2^^.) + l.o.t. 

< ciiDV(t)iu2(^^^(,„))pi(Diw^(t)c)iu2(K^) + \\dtu\t)aLim+io.t. 

< C{e)\\D\u\t)0\\,.^^.) + ||9,n'^(t)C|U.(^,,(.„)) + l.o.t. 

For e small enough, we can absorb the ||^2(-]fj2) term on the right hand 

side of the equation into the left to get 

II^'^'IIl|(b,(.o)) < C7||9i7/(t)ciL2(K2) + /.oi. 

We can now take in time on the time interval [to, to + SR"^] and use (i) to get 

\\D^u''\\L^{[ta,to+SR^],Ll{Bit(xo))) < C{e). 



m 



L'-(iy^''^) For e > 0, 35 and r G (to, to + SR"^] such that for E^[u^{t); B^jiixo)] < e 



and r > 1 , 

sup ||/^^«'^|U^Mo+«^],L^JB«(.o)) < Ci\\u\r)\\^.,EV])- (3.12) 

h 

Pf: To get the bounds that we want for any r > 1, we cannot hope to start our 
time interval at to since all we are assuming is that u^{to) G Hj^. However, from 
(ii), we can find 6o such that 

\\D'^U^\\LmtQ,to+SoR^],Ll{B2R{xo)) < C. 

For any ti G (to, to + ^oR^), then, there exists r G (to,ti) such that 

l2„,h/ M|2 - ^ 



We can now start at t = r with initial data. 

For any choice of (, Uj solving ()3.2p . u'jQ solves the equation 

dtiu^Q) - A'iu^Q (3.13) 
= {l(C,,+i + 2Q, + a-i)(A?z^,) - 2(D+,0 + D_,Q)iD+.u'; + D,,u^) 

Now we choose C to be the cutoff function, 1 in a ball of radius R supported in 



-{ul^.+2ul+ul_,)A\,] 
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B2r{xo). By Duhamel's, we have 

< \\^\t){D\u\r)Cm,. + 1^ ^^_\y/J \D\\s)\'C\\LtiR^)ds + 
+ \\\D'C\\D'u'^\\\Lt(R^) + \\\u'^\\A^C\\\LUR^)ds 
Using Proposition 12. II for s = 2,p = 4, g = 2, we get 

\\D\u\t)0\\Lt 

< mt-T){D\u'^{T)0)\U.+ \,J \\D'u\s)\\h.^^..ds + 



C{R,E'[r]) 



1 



< \\<l>\t-r)iD\u\r)0)\\Lt + C j _ ^y/, \\D\u\s)0\\,. 

■\\D\'\\LliB,,i.,))ds + C{R,EV]) 

Now taking U' in time of both sides, we know by Young's Inequahty that we can 
take ||Z)^u^||2,2(52^(2;o)) -^^ \\D^{u^{s)Q 11^4 in L^; the first term is bounded 
by C(e) by (ii) and the second is exactly the quantity that we are trying to 
bound. This entire product can be absorbed into the left hand side. The linear 
term \\^^{t){D^{u'^{T)())\\^4 is in U since is Hi aX t = r. Therefore, 



L'^{Hl) For e > 0, 3(5 and r G (to, h + ^R^) such that for s < oo. 



E\u\ty, B,n{xo)] < 6 ^ sup \\D'u%LliB^^,,)) < CiWu^T)^^, E^lf^]) 

T<t<tQ+5R? 

(3.14) 

Pf: As in (iii), we can use the L'^{Hl) bounds to find r such that u^{t) G H^. 
For second derivative estimates on u'^C, we need to differentiate ()3.13|) and use 
the estimates (j2.1Up by taking p = 4,q = p' = 4/3. Using Observation 1 and 2 
on and the smoothness of z/, 

Di(A>,) < \D'x';\ + ix'^wdW < \D\i\' + \dW\dW 

To avoid cluttering the basic idea with overwhelming technical detail, we highlight 
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the crucial term: 



\D\u'^m\\Liim 
< 



^ \\CD\X\sMs))\.r.^^.,ds + C{R,E'^[f]) 



< 



h 

\\D\\s) I \D\'^C) Il4/3(M.)^4 + C{R, E^[f]) 

1 

.hi^\/-\\\ , n I rii nl„,/i/'„Ml lllnl„,?i|2/ 



< C\\D\u\r)0\\,.^^^.^ + C j^^jD^^^^^ 
■\\D\u\s)C)\\Limds + C{R,E\f]) 

Now we can take L°° in time of this expression by taking \\D'^{u^[s)Q\\i2^(^^2-^ in 
L°° and ||-D^'u'^(s)||j;^4(^2^(^.g^j in U for r > 4. From this we can conclude L'^iHl) 
bounds in the parabolic cylinder Br{xq) x [r, + 5R^]. 

□ 

Proof. { Theorem \3.^ The construction of a solution to ()3.ip begins, as before, at the 
level of the discrete approximations. Existence for any fixed h is not an issue. Moreover, 
since our initial data is of finite energy, we know from the energy inequality that any 
solution remains of finite energy. However, this energy may concentrate which in turn 
invalidates an immediate argument for the convergence of the interpolants. As seen in 
the case of small energy initial data, the convergence argument depended on having 
control on the second derivatives, or uniform control on the gradient. What Lemma 
13. 31 provides is local control of these quantities. The key extra step that is needed is an 
estimate of the size of the energy concentration set. From this we can conclude that 
the limit of the interpolants can be extended across the concentration sets in such 
a way that it is a partially regular solution of the harmonic map heat flow problem. 

Step 1: Attaining Bounds 

This is the content of Lemma 13.31 

Step 2: Estimating the Size of the Concentration Set 

By invariance of the equation l3.2l bv time translation, it suffices to estimate the concen- 
tration set on the interval [0, 1]. Fix eo and 6 of Theorem 13.21 Since 6 is independent 
of h, it makes sense, for fixed Rq, to partition [0, 1] into intervals Jq, /i, . . . , I2/5RI of 
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length 5RI/2. Denote Tj to be the base time of ly Define /? : [0, 1] ^ {0, . . . , 2/5RI}, 
assigning to every t in [0, 1] the index j such that t & Ij. For arbitrary R and any point 
(xo,to) G X (0, 1] denote by Pji{xo,to) the parabohc cyhnder 

6R^ SR^ 

PRixo,to) = {{x,t)\x e S/j(xo),t G [to ^,to + -1^]} 

Define the singular set 

S = n{(a;,t) G X M+ : limini E^[u^ , PR{x,t)] > eo} 

Then {xk, tk) G S implies, by Lemma that for every R < Ro and every h sufficiently 
small, 

EVitk-^-^),B2Rixk)]>eo/2. 
For any {xk,tk), associate Rk < Rq satisfying 

Then 

E^[u^{Tf3(^t^)), B2R^^(xk)] > eo/2. 

Now these concentration balls have been 'snapped to grid' associated to Rq. At any 
time slice, a direct consequence of the energy bound is that there is a finite number 
of balls of radius 2R in which the energy exceeds eo/2, independent of R. Specifically, 
there can be no more than E^[f^] / (2eQ) such balls at any time slice for any R. The 
balls B2R^ at time Tj cover the bases of the concentration cylinders in Ij in the sense 
of Vitali, so by Vitali's covering lemma, we can take these balls to be disjoint. Now 
adding up the contribution of all the concentration cylinders over the 2/ (SR^) intervals, 
we have 

^kVoliCk) < IsRlJ^kRl < l^^o^j^SRl^^ 

Thus 
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We can conclude that the concentration set S has locally finite 2-dimensional Hausdorff 
measure with respect to the parabolic metric. 

Step 3: Convergence of the Interpolants 

For (xo,to) G S"^, there exists R> such that for some unbounded sequence h, 

E''[u^Pn{xo,to)]<eo. 



27 



Equivalently, 

E^[u\Tt),B^n{x,)]<t,/2, 

so by Lemma IS21 ||-D^'u'^||^2 < C in a uniform neighborhood Q of (xo,to)- Associate 
with each sequence the interpolant p^. From Lemma f2. 31 p'^ — > p in H^{Q). By 
the argument given in the global small energy case we have, for ( G C^{Q), that 

/ dtp-C + Vp-VC-Vp- Vv{p))v{p) ■ C = 0. 

Jq 

Off S, then, p solves ()3.H) . 

To show that p solves the equation on R_|_ x M^, fix a compact set Q CC x M?. 
As in the covering argument above, we can, for fixed i?, partition up the time interval 
occupied by Q into intervals of length 5B?. Without loss of generality, we take the total 
time interval to be of unit length. Given i? > 0, let {Qi^i = PR{xi,l,ti)} be a covering 
of E n Q by parabolic cylinders. Denote by A*"; the number of parabolic cylinders in 
each time interval //. By the energy inequality, A^; is bounded by some finite number 
Aq independent of R, so the number of such cylinders is less than Nq/{5R^). 

Consider smooth cut-off functions < r^i, < 1 identically in -Bi(O) and 1 outside 
of -82(0). Define the cutoff function 

Esr-^ X — Xii I — ti 

1=1 i=l 

Let (j) G C^iQ). Then setting C = 

= / {dtp ■ + Vj9 • V0 - Vp ■ Wiy{p))u{p) ■ (l)}r]R + L 
Jq 

where L can be bounded by 

\L\ < 

< 



< 



< 
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\vpm\vvR\ 



"^Q 1=1 i=l 

~1 R JQ,, 



1=1 i=l 

1 = 1 i=l ^ -^Qr.l JQ^,l 



1 = 1 1 = 1 -^Qi,!- 

Ni 1/&R^ No l/Si^^ 

i=l 1=1 i=l 1=1 

No „ 

i=l 1=1 "^Qi.i 

Taking — > 0, \L\ ^ by absolute continuity so p weakly solves ()3.ip in M+ x M^. □ 

4 LLG 

We proceed to prove our main result, Theorem 11.11 The construction of a solution to 
this problem will follow the same path as that for the harmonic heat flow. The semi- 
discrete approximation for LLG is given by p.2|) . The methods of the previous section 
heavily depended on estimates for the linear equation. One additional ingredient is 
needed for the LLG: a reformulation of ()2.2|1 that yields a system that is linear in its 
highest order term. In the continuous setting, this reformulation is achieved by showing 
that the coordinates of derivatives of u solving ()1.1|) for a = onto a moving frame 
satisfy a system whose highest order term is a linear Schrodinger term. The concept 
of frames is reviewed in Appendix B. For our (discrete) setting, we achieve such a 
formulation by replacing derivatives of u with projections onto the tangent space of 
difference derivatives. 

4.1 The Linearized Discrete System 
Notation 

Difference derivative lose their tangency, so we separate out the tangent and the or- 
thogonal portions. Denote 

where 'V^^^vJ^ G T^^Af and w^^,w^^ ± ^"j-^ • "^^^ tangent portion can be further 
decomposed as follows. We can assume that M admits a global frame {e^, = z/ A e^}. 
Justification for this is provided In Appendix B. Then e* = e*(u^). Introduce the 
complex notation 

i:TM ^ TAf, ie^ = e^, ie^ = -ie\ 
Letting e = and q^'' = q^^'^ + ^g/'''^, qj^ = qj'^'^ + ^Qj^'"^ ■, 
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where A^denotes the orthogonal projection onto .J\f. The orthogonal portion can 
be written as 

jV is a smooth hypersurface so the unit normal is a smooth function. Denote by 

Cz. = ||Vz/||l<x>. 

Since is smooth, e is smooth so we denote Ce = ||Ve||i:oo. 

Observations 

(1) \v+''\ = \q^''\ and \vj''\ = \qj''\ 

(2) 0+^= = D+ku'j ■ Vj = 0(/i|L>+fciiJP). Similarly, af = 0{h\D_ku'}\^). 

(3) Aj = OdD^Mjn. 

(4) For suitable e, and — i^^P < 

This follows from the following manipulation: 

|D+.<r<2(i^;r + Kr) = (|t.;r + |a;«/) 

< 2(i^;|2 + c||<+,-«jnD+,«jr), 

where we have used (1) in the last line. Using a similar argument, for suitable e, 

(5) Since uj = y{u^) 

(6) e is smooth and Cj = e{uj) so 

dtej = V„e9tii^^ = 0{\dtu'^^\), D^ej = C'dL'fcii,^!) 

Lemma 4.1. For e > 0, choose h such that 

l4+i-«JP<e, |iiJ-4_iP<e. 
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Then qj = q^'^ or q - ^ satisfies 

dtq, = iA\j + A^,- + F\qi, D'qi), 

where F^{q, D^q) = 0{\q\^, \q\\D^q\) and I indexes over nearest neighbors of j . 

Proof. Let Uj = Uj, qj = q^^ and Vj = v^^ . To get an equation in terms of g^, take 
-D+fc of ()2.2|) and then take the inner product of the resulting expression with Cj. Let 

/ = D+k{dtU,) ■ e, = dtiD+kUj) ■ e, =11 + III + IV. 

where 

II = A D+kA^u']) ■ ej 

III = {D^,u, A Dl,u^^) ■ e,- 

IV =D+fc(A\^ + Ajz/,)-e, 



= dtivj + OjUj) ■ Cj 

= [dtiqjCj) + ajdtVj) ■ Cj 

= dtqj + qjdtCj ■ Cj + ajc^dtUj ■ Cj 

= dtqj + {ceqj + c^aj)D+kUj ■ '^j)dtUj ■ cj 

= dtqj + (c^qj + c^aj)D^kUj ■ A A^Uj ■ ej + A^Uj ■ ej) 

Each term except dtqj belongs in the nonlinear term F{qj, D^qj). We need to 
show that the magnitude of each term is of the order claimed. 

• {uj A A^Uj) ■ Cj = {vj A D_iVj) ■ Cj + {uj A D_iWj) ■ Cj 

\uj A D_iVj) -Cjl = \uj A {D^iqj + qj^^iD^iCj) ■ ej\ 

= \D-iqj{,Vj A Cj) ■ Cj + Ceqj,-i{vj A vj') ■ Cj \ 

= liD^iqj + Ceqj^-i{i^j A v~') ■ ej \ 

= Oi\D\l\q,\') 



A D_iWj) ■ Cjl = \uj A D^i{ajUj) ■ ej\ 

= \cuaj--i{uj A D^iUj) ■ ej 

= 0{h\D.iUj\^\D^iUj\ 

= 0{{uj^-Uj,^.i)\D^kUj\'^ 

= Oi\q,\') 
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• Similarly, 

\^^Uj ■ ej\ = \D_i{vj + Wj) ■ Cjl 
= 0{\D\\,\qi\'). 

Multiplying both terms by \{ceqj + c,^%)| = the contribution that / has 

to F is 0{\qi\Mqi\\DW). 



II = (uj A D^k^^Uj) ■ Cj 
= {uj A A''{vj + Wj)) ■ Cj 
= (i^.- A AS)- e,- + A AS)- e,- 

• The last term belongs to the nonlinear term F(qj, D^qj). 

I {uj A A'^^;,-) -ejl = i I {uj A - D^i)wj) ■ Cj \ 

= \ I i^j A aj,+iv+' - Vj A aj.^ivj') ■ Cj \ 

= o(^Mi^m>;M^sriv}) 

= 0{\qi\') 

• To handle the term {uj A A'^Vj) • Cj, consider first A'*Vj. 

AS,. = A'^(g,e,) 

= j^{D+i{qjej) - D_i{qjej)) 

= A'^g^e,- + ^{qj.+iD+iej - qj.-iD^iCj) 

The first term yields the linear Schrodinger term 

{uj A A%ej) ■ ej = A ej) ■ ej = iA'^g,-. 
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The magnitude of the second term can be estimated as follows: 
K^j ^ JliQji+iD+iCj - qj.^iD^iCj)) ■ ej\ 

= I \ce{qji+i {^j A vp) ■ ej - q^.^i {vj A vj') • e^} | 

= I \ce{qn+iiq^'" - qj^-iigj") I 
= \ce{qj,+iDo,qj + qj"Do,qj)\ 
= 0{\qi\\DW). 



Ill = D+kJ^jADl^Uj 

1 

= Cr,[D+kUj A -D+kUj+l) 

= ^ ( + ) A + W++ J ) 

= c-\{^J+i A <+i + ^7+1 A v+^^) 

= 7+1 A + w7+i A ) 

= c^^((g7+iej) A {a^+^Vj+i) + (aJ+ii/^+i) A {qj+^ej)) 

= c^^(g7+ia++i(-iej+i) + aj+ig++i(iej+i)) 

This term, which contributes to F, has the correct order since e^+i = Cj + 
and _^ 

|^(g7+ia++i(-iej+i) + a7+ig++i(iej-+i)) • e^'+i] = |). 

= D+k{^^Uj + Aji/j) • ej = A''t;j- • e^- + i^^Wj ■ ej + Aj+iD+fei/j • e^- 
• A'^Vj- • ej = A'^g^- + ^qj^+lD+iej - qj.^iD^iej ■ ej 

\j^iqji+iD+iej - qj,^iD_iej) ■ ej\ = \Ce^{qj.+iV+ - qj^-ivj) ■ ej\ 

= 0{\qi\\DW). 
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-{aj^+iD+iUj - aj^-iD.iUj) ■ ej = C -{aj^+iD+iUj - aj^^iD^iUj) ■ ej 



0{\qi\\D'qi\). 



Aj+iD+fci/j • ej = 0{\qi\ 



□ 



4.2 Proof of Main Result 

Equipped with Lemma l4.ll we can follow the steps illustrated for the harmonic map 
heat flow to prove our main result. 

Lemma 4.2. There is an eo > such that if E^[u^(to); Bi^^xq)] < eo then 36, inde- 
pendent of h such that 

sup l|DV(t)|U2(B_(,„))<C 

to<t<to+5Tf 

Proof. From ()2.3j) . the uniform bounds on energy and ||c?tM'*||i2 are immediate. Multi- 
plying ()2.3j) by dtu'^, summing over j and integrating in t, we have 

EVit)] + \\dtuXHLl)<EV]- 

We will achieve L°°{H^) bounds in the incremental way that was illustrated for the 
harmonic map heat flow: 

(i) A local energy inequality follows exactly as in the proof of 13. 81 using ()2.H|) 



L^(if^) Also using ()2.3|) . nothing changes in the proof of this estimate except 
that dtUj has to be replaced by dtu'j — Uj A dtu^. 



m 



The statement that we will show is: 



For e > 0, 35 and re (to, + ^R^] such that for r > 1 and solving ()2.2|) 

E\u\t)-B,R{X,)] < 6 ^ sup \\D\^\\Lr^^,^,+8R%LiiB,i.,)) < c( || u^r) || , 

h 

Pf: We can use the L'^{Hl) bounds to show that for any ti E (to,^o + ^R^), 
there exists r G (tQ,ti) such that 
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We can now start at t = r with initial data. Since the intention is to use hnear 
estimates on the fundamental solution, we will resort to the linearized discrete 
system given by Lemma HUl Since E'^[u'^{t)] B/^^^xq)] < e, we have in particular 
that 

SO the condition of the lemma is satisfied. Choose to be the spatial cutoff 
function, 1 in a ball of radius R supported in B2r{xq). Then for qj = q^^ or q~^^ 
qjQ solves the equation 

dt{q,Q~{l + z)A'^{q,Q) 
= (1 + + 20 + Q.~i)F\qi, D\) - 

2(D+,0 + D^^Q{D+^q, + D.^q,) - + + ql_,)A%}, 

Applying Duhamel. 

lkWC)llL4(R2) 

k 

< \\U\t - r)(g(r)C)|U4(M.) + _\^,/, {\\M')\' + 

+ lll9l|A'Clll4/3(^.)}rfs 

By definition, g^^'" = Dj^kU^ ■ e^, so 

D^,q^^''' = 0{\D'u%\DW^). 

The \D^Uj\^ term can be seen as the error that appears from using q to estimate 
D^u^; if it were not there, we would immediately have that \D^q\ = 0{\D'^u^\) 
and then exploit the L'^{L\) bounds on to attain bounds on q( from the above 
integral inequality. Instead, we keep in mind that the use of q was purely to 
make use of the linear estimates; it is D^u^( that we need to bound. So long as 
q satisfactorily controls D^u'^, we can achieve the desired bounds. This control 
is a consequence of Observation 4; namely. 
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< lk(*)C)IUt(M^) 

+ ||(|i?V(.)| + |i5V(.)r)5(.)C|L4/3(^.) 



< 



h 

t 



+ |pV(.)|U.(^,,)(,„)||pV(.)CIU4(M^) 



Now taking U in time of both sides, we know by Young's Inequality that we can 
take \\D^u^\\lI{B2r{xo)) (this term is bounded) and \\D^{u^{s)C)\\]^4 in L*". 

Integrating on a time interval of SR"^ and using the LF'^H'I) bounds, the second 
term is easily bounded. Therefore, 



L°^{Hl) Pf: Since 



and wc have from (ill) that has uniform UiW^^'^) bounds in a reduced parabolic 



cylinder, it suffices to attain uniform, local L°°(FF^'°°) bounds on qC^. 

\\D\{m\\L^m 

< \\D'q{T)C\\L^iR^) + (^_^^)3/4 {llg(^)ll^^(M^)lllg(^)l'Cll4(M^) 

+ IP'(k(«)rC)ILv3L + \\i\D'qis)\ + \q\)\D\\lys,^,^} + 
+ C{R, E^[f'])ds 
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< 




+ C{R, E''[f''])ds 



< 




Taking L°° in time of both sides, we have by Young's Inequahty that we can take 
\\'1{^)\\lI{B2r{xo) -^^ > 4 to achieve the bound 



Proof. ( Theorem li.ij) 

Step 1: Attaining Bounds 

This is the content of Lemma 14 .21 

Step 2: Estimating the size of the concentration set 

There is no change in this argument from the case of the harmonic map heat flow. 
Step 3: Convergence of the Interpolants 

The argument has been illustrated in the case of the harmonic map heat flow. The 
only work that has to be done to accommodate for the presence of the Schrodinger 
term is to show that on a compact set Q for which we have uniform L'^{Hl) bounds, 
that the interpolant ph converges. For (f) G C^{Q), we need to show that 



(9t/-0+V/-V0-V/-Vz/(/)0+^{(z/(/)A9i/)-9i0-((9iZ/(/)A9i/)-0} = 0{h) 



Expanding the dtp'^ ■ (j) term, collecting terms in dtu^ and using the L'^{L\) bound on 
dtU^, the only term that is not already 0{h) is dtUj ■ (p which contains the following 
terms: 



It was shown in the harmonic map heat flow case that the terms involving the 



□ 
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damping term have the correct cancellation. It remains to show that 

/ (4 A A\^) ■ + A 9,/) ■ 9,0 - (9,z/(/) A 9,/) ■ 0} 

= /" AD_,^iJ)-Z}+iZ/, AD+i?zJ)-0 + 

(z/(/) A 9,p^) ■ 9,0 - (9,J.(p^) A S^p'^) ■ 0} 

The uniform L°°{Hl) bounds resolves most of the terms. In addition, the following 
manipulation 

/ V D+,{u, A • + (z/(/) A 9,/) ■ 9,0 

JQ i 

combined with these bounds permits us to conclude that ph converges in Q. The 
argument to show that p solves the equation on IR+ x proceeds with no change from 
the harmonic map heat flow case. □ 

A consequence of the proof for Theorem 11.11 is the special case for small energy 
initial data. 

Corollary 4.3. There is a constant eo such that for f G H^iM."^,//), E[f] < cq, there 
exists a smooth global solution to 
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Appendix A 

Let n C and C e C^{n). Then 

|||/rCIU-(M^) < C7(s)||/|U.(n)||V(/C)|U.(M^); '-^ = 1 + 1 (A.l) 

U (J 

Proof. We first take the easy case q < 2. Then by Holder and Sobolev embedding, 

lll/ll/CII|L-(f^) < C\\f\\L,^n)\\fC\\Lrm < C||/|U.(f,)||V(/C)||L.(E2) 
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the last inequality holding iff ^ — | > iff < 2. p, g and s are related as follows: 

1 _ 1 1 _ 1 1 1 
p 2s r 2s q 2 

This is the inequality that we want. 

Now consider the case q > 2. The proof follows in the same manner as that for 
the normal Gagliardo-Sobolev-Nirenberg inequalities, with the help of the following 
identity: 

We can write 

\dl{f'ix,,X2)C{x,,X2))\ 

< 2s r \d,{f{y,,X2)ayuX2))f'-\yuX2)C-\yuX2)\dy,+ 

J — oo 



X-1 



\f\yuX2)d,{C{yuX2))\dy, 



and a similar expression holding xi fixed and integrating on the interval {— 00^X2) ■ 
Multiplying these two expressions together, and integrating, we have 

/•oo /'OO 
/ \f^{x,,X2)C{x,,X2)\'dx,dX2 
- 00 ./ — 00 

< c{s){{r r \v{fo\\p^-'\\r'\dx,dx2f+ 

J —00 J —00 

/CO /*OC /*CX) POO 

/ \V{fO\\f'-'\\C-'\dx,dx2){ / / V\C\\ndxidx2) + 
-00 «/— 00 J —00 J —00 

-OD POO 



/OO POO 
/ \v{c)\\ndxidx2n 

< c(.)||v(/c)||L||/||L|||/rcill^.^ 



□ 



Appendix B 

Let A/" be a compact smooth manifold. Given p G A/", there exists a neighborhood U of 
p and a local coordinate system {y^} on U . An immediate basis for the tangent space 
at any point in U is given by the coordinate vector fields Xi = However, it is 
sometimes not convenient to use vectors that come from a coordinate system. For any 
n linearly independent vectors in TyH 1 we denote the value of the field in p by 

X(p) = {Xi(p),...,X„(p)}. 
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X{p) is a basis for TpAf. A moving frame on A/" is a function F : p E M ^ ^{p)-, where 
F{p) defines a basis for TpM . If F is defined for all p G A/" then F is a global frame for 
M. 

Throughout the process of linearizing the discrete LLG in Section H?T] . we assumed 
the existence of a smooth global frame on M . The problem is that this assumption 
can be thwarted by simple examples, the sphere S"^ being one such case. Fortunately, 
this potentially unfortunate hindrance can be successfully bypassed by at least two 
methods. These methods are standard so we provide here a highly customized sketch 
of these two remedies: (1) embedding Af isometrically as a totally geodesic submanifold 
into a larger space which does admit a global frame and (2) using a puUback frame on 

X M by exploiting the smoothness of the map u. A good reference for this is j3]. 

(1) Embedding into something larger 

One way of introducing a global frame is to consider the equations on a tubular 
neighborhood t/ of A/" in M'^ which does admit such a frame. In this case y & U 
can be represented by coordinates {u{y) , p{y)) . For u G Af, z/(m) can be extended 
to U by z/(?/) = v{u{y)). Then for small p, 

y = u + piy{y). 

The equation for y is given by 

dty = u{y) AAy + au{y) A {u{y) A Ay). 

All derivatives of u can be computed using a global frame on M?. However, this 
procedure has the drawback of losing all the geometric structure of the problem 
and consequently all the estimates and cancellations these structures imply. 

We can view the LLG as a specific example of a Schrddinger map which is a map 
M : M"^ X M (A/", g, J), where A/" is a Kahler manifold with a complex structure 
J = ux and a metric g giving rise to covariant differentiation D, satisfying 

dtu = J{u)Ddu + aDdu. 

What we really need is to embed M isometrically into a larger manifold M which 
itself has a complex structure J where J\j^ = J and in such a way that the 
equation in Af 

dtU = J{u)Ddu + aDdu 

has Af as an invariant subspace. For our problem this can be achieved in the 
following manner. Figure ^ attempts to illustrate this embedding. First we 
embed A/" C C (a complex manifold must be of even dimension.) On a 
tubular neighborhood U of A/" in we introduce the metric g given in coordinates 

u = u + pviu) + ^k, u G Af: 

ds'^ = dslf + dp^ + di'^ 
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Figure 1: Embedding A/" C C M"^ 

and the almost complex structure J{u) 

Jy = J{y{u) + ai'{u) + hk) = Jy{u) — hi>{u) + ak. 

The four dimensional manifold U has a global frame by construction and is given 
by a metric giving rising to a covariant derivative D such that 

Z)J = 0. 

The equation 

dty = JDdy + aDdy 

has as an invariant subspace. This last statement is true since for uq G A/", we 
have Po = ^0) and the equation reduces to 

J dtu = J{u)Ddu + aDdu 

\ U{0) = Uq 

(2) Fullback Frame 

Consider a smooth embedded surface A/" C M'^. A pullback frame for a map 
M : n C — A/" is a choice, for each point x E fl, of an orthonormal basis 
{ei(x), e2(x)} for Tu{x)-^ ■ Since M is embedded, there is a well-defined tangent 
plane T^u{x)^ for every point x G An orthonormal basis can always be specified 
at a point so this type of frame always exists. If M admits a global frame, 
then a pullback frame for u automatically exists, given by {ei('u(x)), e2('u(x))}. 
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However, the converse is not true. The pullback frame is a frame on the domain, 
and might be different at u{xi) and u{x2) even when the two values coincide on 
A^; a pullback frame, therefore, might exist in the absence of a global frame. 

The construction of such a frame is most apparent in the case m : / C M A/". 
Let Co be a frame at T„(o)A/'. Then we can solve for e by parallel transport (i.e., 
e satisfying D^e = 0) which gives the following o.d.e.: 

d^e+ {e- {u{u))^)u{x) = , . 

e(0) = eo ^^-^^ 

Figure El illustrates this selection of a frame. 




Figure 2: Constructing {61,62} by parallel transport 



From this construction it is apparent that the frame reflects the geometry of 
u, since e and u have the same regularity. Higher dimensional analogs of this 
construction are standard. 
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